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On The Reflection and Refraction of Hydromagnetic 
Waves in Ionised Gas 


By Tomikazu NAMIKAWA 


Institute of Polytechnics, Osaka City University 
(Read Oct. 31, 1960; Received March 14, 1961) 


Abstract 


We are interested in the manner of propagation of hydromagnetic waves in 
the Earth’s upper atmosphere. The density changes continuously in the atmos- 
phere and the steady magnetic field is a dipole field, but we discuss here the 
laws of reflection and refraction of hydromagnetic waves at a plane boundary 
between two semi-infinite homogeneous media of anisotropic electrical conductivity 
under a uniform steady magnetic field. When the steady magnetic field is perpen- 
dicular to the boundary plane, the laws of reflection and refraction are simple, 
but when the steady magnetic field makes an angle with the plane, they are 
complicated. The incident ordinary (or extraordinary) waves may be reflected 
or refracted as the ordinary or extraordinary waves and propagate anisotropically. 


1. Introduction 


Some of the geomagnetic disturbances have their origin outside the Earth’s atmos- 

phere. The disturbances propagate as hydromagnetic waves through the Earth’s 
- atmosphere suffering phase delay, reflection, refraction and absorption, and finally are 
observed at the Earth’s surface. Piddington (1956, 1959) provided useful information 
about the transmission of geomagnetic disturbances by a study of the separate equation 
of the ion-electron plasma and the neutral atom gas. 

The laws of reflection and refraction of a plane Alfvén wave in an incompressible 
medium of infinite electrical conductivety have been investigated by Walén (1944, 1946), 
Alfvén (1950), Lundquist (1952), Ferraro (1954) and Roberts (1955). The investigation 
in a compressible medium of infinite electrical conductivity was done by Simon (1958) 
and Williams (1960). We shall discuss in this paper the reflection and refraction of 
hydromagnetic waves under the condition of the upper region of the atmosphere 
(Piddington 1969) in which hydromagnetic waves propagate in the ion-electron plasma 
alone. It is difficult to discuss the propagation of hydromagnetic waves in a dipole 
field. Though this paper concerns weak hydromagnetic waves in anisotropic medium 
under a uniform field, it may throw some light on the manner of reflection and re- 


fraction of hydromagnetic waves in the upper atmosphere. 


2. The transmission equations 


The transmission of hydromagnetic waves through a partially ionized gas is dis- 


117) 
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cussed by using the equations of motion of two gases only, the plasma and the neutral 
atom gas (Piddington 1959), neglecting ordinary viscosity due to thermal motions of 
the particles (Piddington 1957) and the effect of gas pressure gradients caused by the 
wave (Piddington 1955, 1959). The hydromagnetic equations may be listed (Piddington 


1954a, 1956, 1959,) 


pPPh= srcoy| OF p x(vx Hy)| +2 2 (singsiny ~ +sinBcosr 7 + cosh, o \y xh, ‘eh 


ot 
61 5) mea / 
ap hes V+ (Hy x7 xh)=0 (2) 
ov’ el 
a —t@-v)=0, (3) 
where 
4ro 


Here o,, a. and o; are the conductivity components (e.m.u.) of the plasma alone. 
h, H,, v, v’ and n= are the perturbed magnetic field, steady magnetic field (whose 


direction cosines are 0, cosy and cos{), velocity of the plasma, velocity of the neutral 
atom gas and the ratio of the mass densities of neutral atomes and ions. +c is the 
period taken to accelerate the neutral atomes to the same velocity as the plasma. 
Though we are dealing with a partially ionized gas, equation (1) refers to a fully 
ionized gas because the effects of neutral atomes are introduced separately by equations 
(2) and (3). This method is essential when we are dealing with waves of frequency 
greater than the neutral atom collision frequency, when the concept of conduitivity 
in the partially ionized gas is not valid. This method can be applied for most geo- 
magnetic disturbances (Piddington 1959). Equations (1), (2) and (3) are linearized. As 
the perturbations v, v’, and A are small compared with ‘V and Hi, these equations are 
generally valid except at several Earth radii where the field may be violently dis- 
turbed. Plane solutions of equations (1)—(3) are sought of the form exp i(wt—k-r), 
where k is the wave-number vector whose components are (0, sin@, cos@). From the 
equations (1)—(3), we get 


(hk? +tew*)h, +icHycosp kev —ickHywsinBsiny (sindv, + costo.) + 20°F Rhy,=0, (4) 
o,COS 


(hk? + tew* hy + ic cosp kwvy—ickHywsinBcosy (sindv,y + cosbv,)—@2cosd cospk?h,=0 , (5) 
a1 


(k?+iew*)h, + teHocosphuv,—icHywcosa(sindvy +cosdv,) +22 Gand cosph hs =05 (6) 
H,Q- oi APR Gi eo 
H,Q~'cosd wvy +k V*(cosBhy—sindcosdsinfsinyh,)=0, (Cis) 
and 
H,Q-'wv,—kV? (singsiny costh,-+%B0081,,) <9, (9) 
cosé 
where 


cos¢=sin@ sin cosy +cosdcosf , (10 ) 
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_ 4ra® 
e=— Ci.) 
and 
Q= (1+7(1+iwr)}-!. (12) 
3. Dispersion equation 
From the equations (4)—(9), we have 
‘ , 
a =QVS*+—+, (13) 
R € 
where 
1 . Sneed 
S=={1+cos" + — cos)? + 40-608 a 
F ¢ fa cos?)?+ og: (14) 
and 
€0,V*_Q2 
oes (15) 


The term = in (13) may be neglected (Piddington 1959). Because in a magnetic 


field of one gauss 2~10* rad sec™t and w< 10 for all geomagnetic disturbances (Pidding- 
ton 1959), the third and fourth terms of equation (14) may be neglected so that S=cos¢ 
for the O wave; S=1-for the E wave. 


4. Reflection and Refraction of plane waves 


We discuss the reflection and refraction of hydromagnetic waves in the upper 
atmospheric region (Piddington 1959) where the disturbances travel as hydromagnetic 
waves in the ion plasma alone. Losses are small in this region. In the lower atmos- 
pheric region, the medium behaves for waves of all periads between 1 sec and 10‘sec 
as a regid conductor and as a dispersive medium. Let-axes O,, O, and O, be chosen 
in the way shown in Figure 1. Let H, make an angle of 6 with O,, and its projection 
on O,, angle y with O,. Denote the angle of incidence, reflection and refraction by 
6,, 0, and @, respectively. The density changes discontinuously at the plane S. 


incident wave 


refracted /wave 


Fig. 1 Illustrating the reflection and refraction of Alfvén waves 
at the plane surface S. 
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In the upper region wr>1+7 and 


pr? 


aS (Piddington 1959). (16) 
We denote Alfvén wave velocity by V,; and V,, in the region I and II in Figure I. 
We must have, as a consequence of the boundary conditions (A, EH and vy+v/ 

are continuous at the boundary plane S.) which must be satisfied at any time and 


at any point of S, 


Wi= Wr = We » (17) 
and 
ksind;=k,sind,=—k,sind, , (18) 


where the suffix 7, 7, and ¢ denote the incident wave, reflected wave and refracted 
wave respectively. 
a) Incident O wave and reflected or refracted O wave. 
_ From (18), we obtain the laws of refection and refraction of the O wave, 


tand,=tané;(1+2tanfcosytand;)" , (19) 
and 
EOE Bh BR ALS | 
tand,=tandd V+ ( 
We get from (19) and (20), 
6,=+90° for 0,=tan7!{—(2tan8cosr)"'} , (21) 
and 
° = fe Vi\e 
6.= +90 — { Ff ee & ~ 
poet for é;=tan y(t re] (tanf§cosy) F. (22) 


ae 1 )cosytangtand,} (20) 
. 

i) In the case, 8=0 or T=> eqns. (19) and (20) become 
2 


0:=0,, (23) 
and 
V, 
tan@, = Vi tané; . (24) 
ii) In the case, y=0 and A=", we have from (19)—(22), 
tané,=tan0(1+2tané,)~', (25) 
i Ww, ; = 
tand,=tand{ y+(y7'—1) tangy} (26) 
0,=+90°, for  0,=—26°35', (27) 
0=+90°, for  6,=45° (p=2), (28) 
1 
and 
6,=—53°08’ (Yul 
= —53°08 (y=+4)- (29) 


From the curve for the hydromagnetic wave velocity, function of altitude, computed 
by W.E. Francis and others (1959), we see V~10°m/sec., 2X 10°m/sec., and 5X 105m/sec 


F 
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at the altitude of 10*km, 2x10%km and 10?km. We take w= in the region higher 
1 


than 2x10°km and mst in the lower region than 2x10°km when we discuss the 


refraction of hydromagnetic waves. 


Table I. Angles of Reflection and Refraction of O waves when f=0 or r=F- 


0; 90°1...75°.-)..60°% | AS oe Sh P95?) O° PH 158 —S0" ie a5> | 60 | — 75°" —908 
6, 20° 73° OOS | 45") |) 15°-490° |) 2508" — 80° —60° | —75° | —90° 
V 
0:( 73 =2) 90° "amauta wi: al inated woidaed 0° |—28°11'|—49°06'Inadmi-|—73°54’| —82°23'| —90° 
? ssible 
a felg 8°13, 3°50’ 0° |— 3°50’\— 8°13’ —23°25'|—43°01'; —90° 
: | 


l 


Table Il. Angles of Reflection and Refraction of O waves when b=> and 7=0. 


0; 90° 75° 60° 45° 30° 15° 
6, 26°33’ 23°48’ 21°13’ 18°26/ 15" 9°54’ 


—78°04' |—90°, +90°) 69°54’ 36°12’ 


0 ( 4 =2) —63°26' | —69°54’ 


04( a =5) 18°26’ 13°48’ 10°40’ 8°08’ 6°45! alia 
1 
0° —15° —30° —45° —60° —75° —90° 
0° —30° iat 35°06’ 30° 26°33! 
0° —22°55’ | —36°12’ | Inadmi- | —51°44’ | —57°38’ | —63°26’ 
ssible 
0° — 4°47’ | —14°17' _ 5O°22! Papas 18°26’ 


Fig. 2 a. Reflection of O wave when f=0 Fig. 2 b. Refraction of O wave when f=0 or 
or y= y= . Incident waves in the A (B) region refract 
oe 2 


into A’ (B’) region. 
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Fig. 3 a. Reflection of O wave when 7=0 and Fig. 3 b. Refraction of O wave when 7=0, 


ca V, d : 
o> Incident waves in the A (B) region f= and ee Incident waves in the A (B, 
1 . 
reflect into A’ (B’) region and vice versa. C) region refract into A’ (B’, C’) region. 


Fig. 3c. Refraction of O wave when ;=0, 
wee: Va- 1 : : 
ji 4 and lear Incident waves in the A (B, 


C) region refract into A’ (B’, C’) region. 


b) Incident E wave and reflected or refracted E wave 


The laws of reflection and refraction of the E wave are obtained from the equation 
(18), 


sind, =sind, , (30) 
and 


sind,= —Pisind, (31) 
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Table III. Angles of Reflection and Refraction of E waves. 
0; 90° Ys 60° 45° | a0” a 15° : 
a 90° woo 60° 45° 30° 15° 
V 
0, V, =2) —90° —31°10! 
V, 
0.( _ =7) —14°29 | —13°59’ | —12°30’ | —10°11' | — 7°11’ | — 3°49’ 
Von, be 
0° —15° —30° —45° —60° —75° —90° 
0° —15° —30° —45° —60° —75° —90° 
OF 31°10 90° 
0° 3°43’ Colle 10°11’ 12°30" 13°59’ 14°29! 


Fig. 4 a. 


Fig. 4 c. 


Incident waves in the A (B) region refract into 
A’ (B’) region. 


Reflection of E wave. 


Fig. 4 b. Refraction of E wave when 7 me 
¥ 


Incident waves in the A (B) region refract into 
A’ (B’) region, while incident waves in other 
region are reflected back. 


Refraction of E wave when 


rhea 
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c) Incident O wave and reflected E wave 
The law of reflection is 


sin0,= sate =tan0;(sinB cosy tan@;+cosf)"!. (32) 
c 


i 


As |sin@,|<1, we get, 


tan~( —cosf )<<tan( cos ), fair: cosf>0, ), (33) 
1+sinfcosy 1—sin fcosr and 1—sin2fcos*y 0 
and 
a.<tan-*( —cosB ) feseee cosp>0 ) F (34) 
2sinfcosy and sinfcosy=—1 


If B=r=4, we get from (33) 


—25°14'< 6;< 54°44’ (35) 


Table IV a. Angles of Reflection of E waves when the Incident waves are the O 


waves and p=r=4- 
6; s4caa’ | 45° =| 30° 15° | 0° | —15° | —25°14’ 
0, 90° 55°56 | 35°26" 18°33’ o° | —25°03’ | —90° 


Fig. 5 a. Reflection of E waves when the 
Incident waves are the O waves and 6=r=4- 


Incident waves in the A (B) region reflect into 
A’ (B’) region. 


y 


d) Incident O wave and refracted E wave 
The law of refraction is 
sing, = —Vusinds __ Vis 


1 . : 
V,cos¢; V, tané,(sin8 cosy tan@;+cosg)-!. (mes 
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From the condition |sind,|<1, we get 


—cosp —cosp 
tan- Vi1<6,m tan r ; 
1 


sinB cosy +— i Sec 
T Vi sinf cosy V 


when cosf>0, 
: i 
and see aN | ey 
V1 
P po —cos 
;>tan- tan| Vi. kT 
v7, ~sinBcosy or @<tan Vt 4 sin BCosy , 
1 Vi 
when cosf#>0 
38 
and sin®gcos'y—T")'s0 a 
Vi 
when cos#>0 
and 6;< tan-1 Vicos8 V (39 ) 
2Vi1 and sinfcosy=——"! 
if B=r=7) we get 
° / ° / Vis _ 
—15°47'<0,< 25°14’, (sil 2h (40) 
1 
°4Q/ 999/ Vir_ 1 
and 6;>—43°19 or 4,< —70°32’, ory (41 ) 
1 
Table IV b. Angles of Refraction of E waves when the Incident waves 
ae 
are the O waves and f=;=—, y=: 
1 
6; 25°14’ id? 0° —15° —15°47’ 
Oo [=> =2) —90° —39°34’ 0° 69°12’ 90° 
1 


Table IV c. Angles of Refraction of E waves when the Incident waves are the O waves and 
of nt 


6; 90° | 75° 60° | 45° 30° 15° 
a,(+=4) 39° | —21°15' | —15°59 | —11°57" | — 8°20’ | — 4°34 
4 
o° |-—15° | —30° . | —43°19” | —70°32" | —75° 


0° 6°43! 20°11’ 90° —90° —53°37/ 
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Fig. 5 b. Refraction of E waves when the Fig. 5 c. Refraction of E waves when the 


» a 
Incident waves are the O waves and p=r=F> Incident waves are the O waves and f=y;=— 


4 > 
V. 1 ; : ; 
Tp, Incident waves in the A (B) region 7 a: Incident waves in the A (B, C) region 
1 1 ‘ 
refract into A’ (B’) region. refract into A’ (B’, C’) region. 


e) Incident E wave and reflected O wave 
The law of reflection is 


tand,=cosfsin6;(1+sinf8 cosy sin@,)~! . ( 42 ) 
If p= r= this becomes 
tan0,= A 2sind; : ( 43 ) 
2+siné; 
f) Incident E wave and refracted O wave 
In this case, we have 
V. ; Vinx ; =4 . 
tan@, =iicospsind,(1 —7pisingcosrsiné;) , (44) 
when B=r=7 we get 
tand, = /2 Musing, (2—THsino,) * (45 ) 
Vi a 


Table V. Angles of Reflection and Refraction of O waves when the Incident waves are the EF 
waves and S=;=~., 


4 
0; 90° ie 60° 45° 30° 15> 
6, 20° Lat 24°44’ 23°08’ be Ara hf 15°47’ Shea 4 
Vir a 
O, lon =2) 90 88°34’ 83°45’ 73°39! 54°44’ 20°17" 
Vu me! oor 
o( V, =5) 125 11°00’ 9°45' 7°48! iageay 2°42f 
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oF =" —30° —45° —60° | 75° —90° 

0° —10°52" | +25°17" | —37°43’ | —47°12" | —52°38" | —54°397 
0° —16°13' | —26°09’ | —30°22’ | —33°17' | —34°48’ | —35°16/ 
0° = 2°82" | — 4930! 9 — 6°38) | =):7°52" |; 8°40’ |. — 8°56" 


Fig. 6 a. Reflection of O waves when the Fig. 6 b. Refraction of O waves when the 


Incident waves are the E waves and p=r=4- Incident waves are the E waves and f=;=~, 
Incident waves in the A (B) region refict into 7 2 Incident waves in the A (B) region 
1 
A’ (B’) region. refract into A’ (B’) region. 


Fig. 6 c. Refraction of O waves when the Incident 


tus me Vuhip ly 
waves are the © waves and f=7= rie al 
Incident waves in the A (B) region refract into A’ 


(B’) region. 


5. Concluding Remarks 


We can see from Fig. 3—Fig. 6 the complexity of the laws of reflection and 
refraction of hydromagnetic waves when the steady magnetic field H, makes an angle 
with the boundary plane. When H, is perpendicular to the boundary plane, the 
incident O (E) wave may be reflected or refracted as O (E) wave, but when A, is 
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inclined to the boundary plane the incident O (£) wave may be reflected or refracted 
as O and E waves. In the upper atmosphere, the steady field Hj is a dipole field and 
density changes continuously. Under these conditions, it is difficult to discuss the 
propagation of hydromagnetic waves. This paper discusses it in a uniform field and 
concerns only with the laws of reflection and refraction. A complete determination 
of the reflected and refracted waves requires the calculation of the amplitudes. Further 
investibation on this problem may be useful to some extent for the interpretation of 
the observational data of geomagnetic disturbances. 
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Locality of Calm Day Atmospheric Electric 
Potential-Gradient in Kyoto City 


By Akira OKAWATI 


Geophysical Institute, Kyoto University 
(Read May 16, 1960; Received March 17, 1961) 


Abstract 


The simultaneous and continuous observations of atmospheric electric potential 
gradient were made at several stations in Kyoto City which has about one million 
population. The diurnal variations in calm days were observed, and the charac- 
teristic diurnal variation in Kyoto City as well as the locality within some kilo- 
metres were clarified. The diurnal variation at Iwakura 6km apart from the 
centre of the city had quite a different type from the stations near the centre. 
The station 5km. apart from the centre had a diuranal variation partly influenced 
by the centre. And the station 4km. apart from the centre and about 120 metres 
higher than the others had the different feature from others. 


1. Introduction 


It is generally considered that the earth’s surface makes a concentric spherical 
. condenser with the surface of the earth as the negative plate, the upper conducting 
laver or equalizing layer (it is estimated more than 20 km. above the ground but some 
investigaters considered it as the ionospheric E layer about 100 km above the ground), 
as the positive plate and the atmosphere as the dielectrics. Since the atmosphere is 
ionized by the cosmic rays, the atmospheric radioactive substances and the radiation 
from the radioactive substances in the earth’s surface, the atmosphere has finite electrical 
conductivity. The leakage of electric charge through the condenser gives an air-earth 
conduction current. The calm day diurnal variation of atmospheric electric potential 
gradient observed at the ground is mainly due to the variations of air-earth current 
and of the conductivity of air near the ground. 

It became clear by the investigations of C.T.R. Wilson (1920), S.J. Mauchly (1923), 
F.J.W. Whipple (1929), T.W. Wormell (1927) and O.H. Gish and G.R. Wait (1950) that 
the electric potential of the conducting layer is maintained and varied by the thunder- 
storm activity over the globe. So the diurnal variation of atmospheric electric potential 
gradient at the place where there is little variation of air conductivity should have 
the same phase as the electric potential of the conducting layer. 

It was found by the observations during the cruise of “Carnegie” made by J.P. 
Ault and S.J. Mauchly (1926) and the observations at Jungfrau and Zugspitze by H. 
Israél (1951) that over the ocean, at the polar region and on high mountains where there is 
little variation of conductivity, the atmospheric electric potential gradient varies in the 


(129 ) 
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same phase depending upon the universal time anywhere over the globe. 

On the other hand, on land the air conductivity is varied by air-pollution and 
ionizing agent in air, hence the diurnal variation of atmospheric electric potential 
gradient contains large local effect depending on the local time. J.G. Brown (1935) 
investigated the diurnal variation of atmospheric electric potential gradient based on 
the results of observations of twenty-two stations in the’ world. He called the variation 
on universal time the “ unitary variation”, and he considered that the variation over 
the ocean was composed only of the unitary variation, but on the land, the local variation 
superposed on the unitary variation. He tried to separate the variation on land into 
the local variation and the unitary variation. Comparing the local components of each 
stations, he concluded that the local variation of atmospheric electric potential gradient 
is due to the variation of electrical conductivity of air. 

As it is well known, the conductivity of air (A) is expressed as follows; 


A=nek (ig 


where u denotes the density of small ions in air, e the electrical unit charge quantity 
and k the ionic mobility. Small ion is produced by the radiation of cosmic ray, atmos- 
pheric radioactive substances and terrestrial radioactive substances in the surface of 
the ground, and is destroyed by recombination with a small ion of opposite sign giving 
neutral molecules or neutral nucleus and attachment to a nucleus giving a large ion. 


It is expressed simply, 


a =q—an*—BnN ( 2 ) 


where m is the number of small ion, g the rate of ion pair production, a the coefficient 
of recombination, f the coefficient of attachment and N the number of nucleus and the 
large ion of opposite sign. Over the ocean N is ignored, being compared to ~. Then 
the equation becomes in the equilibrium state, 


q—an?=0 or n=,/4 ( 3.) 
a 


On land, when N is much larger than ” the term an? is ignored. Then in the equi- 
librium state the equation (2) becomes, 


—BnN=0 he =— : 
q B fe) n = (4) 
From (1) and (4), we may write, 
1 


The electric potential gradient is expressed as, 


E=-- (6 ) 


Where ¢ is air-earth current. From (5) and (6), we get 


EcctN € By 
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The atmospheric electric potential gradient is roughly proportional to the number 
of nucleus when the air-earth current is constant. 

F.J.W. Whipple (1929) observed the potential gradient as well as the air pollution 
and found that the phase of both variations coincides well with each other and that 
the clear displacement of the phase of the diurnal variation of potential gradient 
appeared after the introduction of Summer Time. This means the human activity in 
large city influences much on the atmospheric electric potential gradient. 

R. Mihleisen (1955) studied the potential gradient and conductivity of air in large 
cities paying attention to the influence of various industrial and artificial equipments 
on the elements of atmospheric electricity. He used portable instruments and measured 
the elements at such places as gas-works, railway stations and nearby streets. He 
found that the smoke, air pollution and exhaust gas of engine are not electrically 
neutral but have charge initially. He also measured the electrical charge produced at 
combustion and found whether the charge is positive or negative depends on the kinds 
of substances burnt and the condition of combustion. He concluded that the variation 
of potential gradient could not be explained only by the variation of conductivity of 
air, and that in large cities the positive space charge is much produced and tends to 
accumulate easily. Also the fact that the potential gradient on Sundays is much less 
than in the week days shows clearly the industrial and artificial influences on the 
atmospheric electricity. 

In large cities, even on calm days the atmospheric electric potential gradient varies 
remarkably and has an interesting feature of variation. In order to study the locality 

_ of the potential gradient in Kyoto City and her suburb, the simultaneous and continuous 
observation was carried out near the centre of Kyoto City and in her suburb. 


2. Instruments and Observation Stations 


In the first period of observation between December 1955 and March 1956, a 
mechanical collector was used, and it was connected to a recording milliammeter through 
the direct current amplifier having high input impedance. In the second period of 
observation between January 1959 and May 1959, radioactive polonium and a relay of © 
special type were used as a collector and it was connected to the same amplifier and 
the recorder. The measurement of the absolute value of the potential gradient was 
carried out in a suitable plane as near as possible to the stations to get the factor of 
plane reduction. 

In Fig. 1 the observation net, its topographical environment and the industrial 
and thickly populated area are shown. In the first period, the observation stations 
were selected at A, B, C, D, E and F shown in the map, and in the second period, 
S047.DiwG,) Hy tJ rand K. 
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Fig. 1. Topographic map showing the distribution of observation stations and industrial 
and populated area of hatching. 


3. Results of Observations and Discussion 


From the data obtained in the first period between Dec. 1955 and Mar. 1956 the 
calm days were selected and the hourly values referred to the Japanese Standard Time 
(135°E or 900") were averaged for about three months. The results are shown in 
Fig. 2. Though these observation stations are not far more than 4km from each 
other, their mean values as well as the shapes of diurnal variation curves are 
remarkably different. The values at about 4” are almost all the same, but the nearer 
the centre of the city the stations are, the larger the mean values are and the broader 
the amplitudes of diurnal variations are. It seems to the author that the air is clean 
at all the stations towards 4”, but later in the morning the nearer the centre of the 
city the air is, the more it comes to be contaminated. Though the potential gradient 
begins to increase rapidly at the stations near the centre of the city at 5", at C station 
it begins to increase at 6” and the rate of increasing becomes gradual towards 8", 
then again increases rapidly. The farthest station from the city is D, where the potential 
gradient increases gradually in the morning. Both the C and D stations have obscure 
second maximum in the afternoon. These characteristics at C and D stations are due 
to the fact that they are distant from the centre of the city where the air is contaminat- 
ed much more. 

More precise observations were needed in order to study the locality and the 
structure of the atmospheric electric field in and near Kyoto City, and the observations 
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were carried out again in the second period between Jan. 1959 and May 1959. The 
calm days’ hourly mean values were summarized for each month and the results are 
shown in Fig. 3, 4,5, 6 and 7. As it is well known, the potential gradient is larger in 
winter and smaller in summer. The results show that the nearer the centre of the city 
the stations are, not only the larger the mean values and the broader the amplitudes 
of diurnal variations are but also the more remarkable the seasonal variations during 
this period are. The seasonal variation is little at G station 6km apart from the 
centre of Kyoto City, moreover the shape of diurnal variation curve varies little. The 
diurnal curves obtained in the second period have the second maximum less than in the 
first period even near the centre of the city, and they come to have obscure ones 
month after month. 

The diurnal curves at D station show rapid increase towards 6” coinciding with 
the stations near the centre, but the rate of increase becomes gradual towards 7” then 
again rapid. This is quite similar to that of C station in the first period. The potential 
gradient begins to increase rapidly towards 6” at G station coinciding with the stations 
near the centre, but it dose not reach the value twice as much as the value of 5” even 
at maximum not like the stations near the centre. The amplitude of diurnal variation 
at G station is much less than others. Let’s compare the diurnal variations of G and 
D stations. We find they are similar in the course of diurnal variation except the part 
between 8” when D station again begins to increase rapidly and 12”. This suggests 
the electric field has twofold structure at D station; the rapid increase towards 8” is 


O 3 6 9 12 1§ 18 2\ 24 h 


Fig. 2. The diurnal variations of atmospheric electric potential gradient in the first period 
between Dec, 1955 and Mar. 1956. 
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Fig. 3. The diurnal variations of atmospheric electric potential gradient in January in the 
second period between January 1959 and May 1959. 
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Fig. 4. The diurnal variations of atmospheric electric potential gradient in February in the 
second period between January 1959 and May 1959. 
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Mar. 1959 
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Fig. 5. The diurnal variations of atmospheric electric potential gradient in March in the 
second period between January 1959 and May 1959. 
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Fig. 6. The diurnal variations of atmospheric electric potential gradient in April in the 
second period between January 1959 and May 1959. 
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Fig. 7. The diurnal variations of atmospheric electric potential gradient in May in the 
second period between January 1959 and May 1959. 


due.to the influence of the centre of the city where industrial and social activities are 
vigorous and so the air is much contaminated. And if there should not be large 
industry and population in Kyoto City, D station would have the same mode of diurnal 
variation as G station. The contaminated air produced much in the centre of the city 
diffuses around and begins to influence D station towards 8”. The two stepped increase 
observed at C station in the first period would be explained so. It is imagined that 
if there were no industrial and artificial influences in Kyoto City, all the stations would 
have the same diurnal variation as G station. 

At H station the potential gradient begins to increase about one hour later than 
others. This is due to the location of H station which is on a hill about 120 metres 
higher than other stations. The contaminated air having less electric conductivity and 
much positive space charge which is produced on the ground is prevented to flow 
upwards by the inversion layer which often appears early in the morning, and the 
potential gradient on the hill dose not begin to increase till the inversion layer comes 
higher or vanishes away. 

On Feb. 3, 1956 it was very cold in the morning, air temperature was below 0°C 
untill 8”, it increased as much as 7.3°C of maximum at 15” and it drew near 0°C 
towards 24". It was cloudy all day long and there was no wind except in the time 
from 12" and 15" during which the wind had the speed of about 2 m/sec. The potential 


gradient of all the stations had the diurnal curves of single maximum type, and had 


no depression towards 15” not as usual. In this day the concentration of contaminated 
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air did not decrease in the afternoon neither by the convection caused by the heating 
of air from the ground nor the diffusion by wind, and it was accumulated near the 
ground. It is imagined that on any day having little wind and little solar radiation 
all day long, the “depression” so called by J.G. Brown, (1935) do not superpose on the 
diurnal course of potential gradient. On the days when south-westery wind is dominant 
C and D stations have usually large potential gradient. This means the influence of 
the city is conveyed directly by the wind. 


4. Conclusion 


There is remarkable locality of atmospheric electric field in the area within 6 km. 
in and near Kyoto City. 

The stations near the centre of the city have larger mean values, broader amplitudes 
of diurnal variation and seasonal variation of atmospheric electric potential gradient 
than the stations far from the centre. 

Within a certain distance from the central area of city, the potential gradient is 
affected by the central area which have. much industrial and artificial contamination 
in the air. 

On the hill about 120 metres higher than the flat, the rapid increase in potential 
gradient lags behind about one hour. 

It seems to the author that the simultaneous and continuous observations not only 
of the atmospheric potential gradient but also of the conductivity of air, space charge 
density and of air-earth current would afford us more precise knowledge of qualitative 
and quantitative character of the diffusion of contaminated air in the lower atmosphere. 
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Abstract 


On the discussion of the charge on rain drop, there remains a serious 
problem on the electric field below cloud, that is the discrepancy between the 
Altielectrograph results by Simpson and Scrase and Simpson and Robinson and 
the theoretical calculation by Whipple and Scrase on the basis of Wilson’s Space- 
Charge Theory. Taking the life time of small ions into consideration, it was 
found that there is the possibility of rapid increase of the field intensity below 
the level of 100m above the earth’s surface ; the magnitude at that level is about 
ten times as large as that at the earth’s surface. This rapid increase may have 
caused the erroneous estimation of the coefficient in the relation between the 
electric current through the sounding apparatus and the field intensity used by 


Simpson and Scrase, and so the above discrepancy may have arisen. 


1. Introduction 


On the electric field intensities at points in the air below electrically active cloud, 
Whipple and Scrase (1936) had calculated the expected values, according to Wilson’s 
space-charge theory (1925), involving the space charge of small ions produced by point 
discharge. Unfortunately the records of the Altielectrograph, sounding records, by 
Simpson and Scrase (1937) and Simpson and Robinson (1940) showed little variations 
of the field intensity, though expected to show noticiable increase. There arose a 
discrepancy between the results of measurements and the results deduced from the 
theory; the former showed nearly uniform field throughout the region below clouds, 
while the latter demands appreciable increase in magnitude with height, resulting in 
the ratio of the magnitude near the cloud base to that at the earth’s surface to be 
ten or more. This has been a problem in the general theory of the electricity in 
stormy weather condition. 

The authors have investigated the rain electricity, and think this discrepancy serious 
for the explanation of the electrical phenomena in rainy weather. One of the interesting 
and striking phenomena in the rain electricity is the inverse-relation between the rain 
current and the potential gradient when the field at the earth’s surface is sufficiently 
strong to produce appreciable point discharge. Simpson (1949) obtained, arranging 
his continuous observational data, a relationship between them and suggested that the 
relation could be accounted for, if rain drops acquired their charges by the process 
of capturing small ions being produced by point discharge and moving upwards in 


(188 ) 
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the field. Such ion capture being considered as the process proposed by Wilson (1929), 
Simpson pointed out that the theory could not be satisfied unless the field intensities 
at higher levels were several times as large as that at the earth’s surface. All of the 
other studies (Hutchinson and Chalmers (1951), Chalmers (1951), Smith (1955), Ogawa 
(1960) and etc.) on the rain electricity demand also the increase of the field intensity 
with height. It is, therefore, necessary to have more complete knowledge about the 
electrical condition below cloud and to resolve the discrepancy, in the discussion of 
small ions’ effect on rain-drop charge. 

It is considered that there are three possibilities of resolving the discrepancy; (1) 
there are some factors omitted in Wilson’s space-charge theory which would reduce 
the space charge below cloud; (2) the analysis of the Altielectrograph record was based 
on too simple relationship between the field intensity and the point-discharge current 
to have information about the field in free air, though it was planned for rough esti- 
mation of the field intensity; (3) the small ions released by point discharge are not 
sufficient to deform the field; this requires much larger effective seperation of the 
discharging point than has been-accepted. Various attempts were made on the first 
of the above possibilities; the effect of recombination process of small ions, vertical 
air-current or the charge on rain drops (Chalmers (1939, 1944)), but no striking efficient 
factor has been proposed. On the second posibility, there had been, for a long time, 
no exact relationship between the two elements. Recently Chapman (1955, 1958) and 
others improved its investigation and showed details about it. 

This paper is mainly concerned with the first possibility and somewhat extended 
_ to the second. The small ions released into the air by point discharge may attach to 
condensation nuclei and others to decrease the concentration with the life time t and 
lose their mobility, then some of them can not proceed above a level below cloud. 
This process may be an important process for the discrepancy in the electrically steady 


state for some minutes. 


2. Electric Field Profiles Expected from Theoretical Consideration 


In order to provide a quantitative treatment of the phenomena, it is necessary to 
make some assumptions and approximations which are here set forth. The theoretical 
calculations are based on the ideas of wide spreaded cloud and so horizontal equi- 
potential surfaces. The field is deformed by the space-charge effect of small ions with 
one sign which are originated in the air by point discharge near the earth’s surface. 
The electric current flows as a result of the small ions’ motion along the lines of force, 
in the absence of appreciable wind. The electrical condition is, furthermore, assumed 
to be steady, involving no time variation. As for points originating discharge, they 
distribute apparently irregularly, but may be considered as uniform sources of small 
ions except quite close to the earth’s surface. 

(1) Wilson’s Space-Charge Theory 

Before proceeding discussion, it is convenient to comment on early theories. 
Wilson’s space-charge theory is the fundamental theory, and the results of the other 
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modified theories differ little from that of Wilson’s theory. 

Wilson assumed that the small ions with one sign, produced by point discharge, 
are the only carriers of electric current and space charge. They ascend with a constant 
mobility in the strong field. 

Let J be the total current density per unit cross sectional area, 2 the concentration 
of small ions of charge e and mobility w and X the vertical component of the potential 
gradient, being put positive for convenience. On the calculation the MKS rationalized 
unit is here refered to. Then, in steady condition, J is the same for every height x 


and 
T=neoX . Cr 
Applying Gauss’s theorem 
dX 
SS He 2 
“ts ( ) 
thus, we obtain the field profile that 
— 21 
SC dearar S (3) 


where X, is the value that would have for y=0, at the earth’s surface. Eq. (3) is the 
result given by Wilson (1925) and usually called Wilson’s space-charge theory. The 
concrete calculation for Eq. (3) was given first by Whipple and Scrase (1936) who 
obtained a numerical relationship between J and X, on the basis of their observation 
on the current from the point connecting to the earth. The numerical relationship 
between the two elements is expressed as 


I=72(Xe—-M?), (4) 


where L is the effective seperation of the reference point with the estimated value 
of 25m. M is a minimum field intensity below which no discharge occurs. a is a 
numerical constant. For positive current entering the point, M is 780 V/m and a is 
8x10-'*; for negative current, M is 860 V/m and a is 10X107'*. The approximate 
results based on their relation are shown by dotted curves in the figure for some 
values of Xo, giving a quite considerable increase of the field intensity with height. 
The problem mentioned in Section 1 is fundamentally the antagonism between the 
Altielectrograph results and the profiles shown by dotted curves in the figure. 

(2) Effect of Immobilized Ions on Wilson’s Theory 

Now we shall consider the behavior of small ions and their effect on Wilson’s 
space-charge theory. Although some of the factors neglected in Wilson’s theory have 
been examined for their effect, being expected to be important, they could not give 
noticiable different features from those given by Wilson’s theory. The attachment 
process of small ions to large particles being roughly considered too (Chalmers (1951)), 
it has thought to be negligible because of too short travelling time of small ions 
compared with their life time. It may, on the authors’ opinion, be more profitable to 
treat of the attachment process by Eulerean method rather than Lagrangean method. 
Some of the small ions originating in the air near the ground may be immobilized 


lat tit 
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on large ions as the result of attachment to condensation nuclei and others. The large 
ions may be accumulated with time and expected to become unnegligible for their 
space-charge density. 

Consider the concentration n of the small ions in Wilson’s theory, then from Eqs. 
(1) and (3) » can be written as 


div teh 
ta to [Nee nt Xe. (5) 


enw 


Pt 


While in the case that the attachment process is proceeding, the time variation of 
small ion concentration at a height x,, is denoted as 


= — 2 (noX)—*, (6) 
where ¢ is the life time of small ions. At the same time space-charge density o at 
%, varies with the rate of 

OT =~ py neoX). (7) 

The source of small ions is limitted in the spatial point near the earth’s surface. 
Its efficiency varies with the field intensity in the way such as expressed in Eq. (4). 
Eqs. (6) and (7) mean the decreases of current in the region above xy, and of the field 
intensity below z;, under a constant electric activity of the concerned cloud; the current 
and the electric field intensity will decrease throughout the whole region, even if the 
initial state were held as expected in Wilson’s theory. 

To estimate the time scale of their reduction, it is assumed that the concentration 
* of small ions produced by point discharge is given in two ways; for the consideration 
of attachment process, it follows Eq. (5) with the exception that the denomination on 
the right-hand side remains constant, and for the consideration of space charge or 
current composed of small ions, it follows just strictly Eq. (5). These give the virtual 
limitting case in which the reductuon rate is underestimated than in the vigorous 
treatment. The transportation of large ions by means of electric force, air-current 
and turbulent diffusion are rather important in the following discussion than in the 
discussion in Wilson’s theory. Thus on the concentration for attachment process ” 


is assumed to be given by the equation 


ctoray Leelee oe 21 o* 
n= X* = cw man a+Xor, (8) 


where * indicates the initial value which is given by Wilson’s theory, and on the 
concentration for the effect of small ions themselves, Eq. (5) is assumed to hold. 

The increasing rate of total space charge Q in a unit cross sectional column between 
a certain height and the earth’s surface is expressed as 


ag éF, ary 9 
ame | ee (9) 


where P, and P, indicate the space charge sum carried on small ions and large ions, 
respectively, in the column. P, is obtained by integrating me from 0 to h on x, so 
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that 
P=" ne-dy ; (10) 
0 
From Eqs. (5) and (10), we have 
dP, _dXo d Wa FDL } (11) 
dt dt dX ih . ; 


By similar treatment, using Eqs. (6) and (8), we have 


See alco De | ite 
dt ct I* ae RS a te tea 


where the first term on the right-hand side in Eq. (6) vanishes because of the de- 
pendency on Wilson’s theory. Applying Gauss’s theorem, X, and X,, the field intensity 
at y=h, are related by Q which takes here negative, as 


Xo=Xn—e0oQ. (13 ) 
If X, is set as undeformed applied electric field, directly produced by the constant 
electric activity of cloud, then the time variation of X, is written as follows, substituting 
Eqs. (9), (11) and (12) for the corresponding terms in the differential from of Eq. (13), 


dX “| 2I ; j=-24f 2I } 
dt dX, zene x +C coon r T*® a re ht Xo Xo* ; ( TH) 


where C is the integral constant. Accepting Eq. (4) as the function of J, Eq. (14) gives 
the solution 


Ma DEE EMM Cexp(— 2KM_t (15) 
V(@+1)X,?—a2M?2+M (a?+1) 5 
or 
an): a 1 KM t 
ob. |i os baked eee 
(a0 a®+1 atl ( ed =?) npey 
where 
2 _2ah 
¥ ~~ egw? 


1 
K= x my Ver 1)X,**—a®M*—X,* } 


and C) or @ is the integral constant which would be determined by the initial condition. 
When we take the numerical values that o=1.5x10-!m?/V-ses and h=10°m 
ordinary acception, and a=9xX10"', L=25m and M=820 V/m, deduced from the 


relationship of Whipple and Scrase, then a? has the value of 2.17x102. And so Eq 
(16) can be written approximately 


Ay Peed KM ,t 
(5g) = oot ea (7) 


therefore, accept the following as approximate equation instead 
of Eq. (17). 
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PE ame 1 oy ee ae 
ay) =1+7,coth*( X +8), (18) 


Oe Crna eZ 


and in particular when X,>10,000 V/m, ie. X,*SM, 


Yon~ x(t). (18’) 


We can not recognize the state to be steady unless it varies in little percentage 
fluctuation about a certain mean state in the time intervals of minutes. Since the life 
time of small ions, tr, may be considered about 30sec on land stations, Eq. (18) or 
(18’) means that Wilson’s theory can not describe the steady state for the strong electric 
field in which the small ions produced by point discharge appear to be appreciablly 
effective. 

(3) Proposed Electric Field Model 

Under the condition above mentioned, space-charge increase in the form of large- 
ion production by attachment process destroies the state expressed by Wilson’s theory 
and reduces the electric field intensity at the earth’s surface to the critical intensity 
below which point discharge occurs hardly. In order to retain the strong electric field 
in steady state, it is necessary to introduce another kind of current which can offset 
the space-charge increase. The ordinary conduction current, which forms the current 
when no point discharge occurs, loomes up as a possible one. It is expected to have 
an important role at higher levels near cloud, though negligible small compared with 
point-discharge current density at lower levels near the earth’s surface. 

Now we shall try to obtain the steady state in concrete features in strong electric 
_field. The assumptions and approximations are taken as same as those in previous 
sections. 

Consider the electrical elements in a unit cross sectional column. Total current 
I, is expressed as 

L=i’X +neoX , (19) 


where J is used as the notation for the ordinary atmospheric conductivity, that would 
be in existence without point-discharge current. The first term on the right-hand side 
in Eq. (19) expresses ordinary conduction current and the second term does the con- 
duction current composed of the negative small ions’ motion, originating from point 
discharge. For the value of 4 without detailed data in the concerned region and cloud, 
it is assumed to be constant with height. The time variation of concentration of 
negative small ions produced by point discharge is expressed at any level in the region 


except the boundary that 


Sf eee eX) (20) 
In steady state 0/,/0x and 0n/Ot are zero. Thus, Eqs. (19) and (20) can be written 


in differential forms, respectively, as 


dX de ut 21 
Vig +O dg MX) =9: (21) 


144 S. SAGA and T. OGAWA 


—of (nX)—2=0. ==’ 


Combining Eqs. (21) and (22) with (19) 


CD Ge Ga ee 23 
Maul wakp ees (23) 
The integration of Eq. (23) is 
[sda Aix 
be eee | -4{4- iy Bee i]. (24) 
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Eq. (24) is calculated for the values at Kew, the station at which the discrepancy 
arose, adopting the same values as those in Section 2 and assuming 4=2x107'* and 
t=30sec. It is noticiable feature that the mean values of the atmospheric electric 
elements at Kew are appreciablly different from those at other usual land stations, for 
example the mean potential gradient in fair weather at Kew is about 400 V/m, the 
value 4 times as large as those at usual land stations. This suggests the necessity to 
consider whether the assumed value of 4 and r are suitable for those at Kew or not. 
The time variation of concentration of positive or negative small ions in fair weather 
follows the equation, 


ane ae 
Tada men, (25 ) 
or 
ang” (25') 


where q is the production rate of small ions, a recombination coefficient, 7 attachment 
coefficient and N concentration of large nuclei. On land stations the second term on 
the right-hand side in Eq. (25) can be neglected. Thus 


d. 
p= I— IN ( 26 ) 


and 

1 oN 

lay (27) 
are approximately hold. Since the atmospheric conductivity may be given by 


A=2neu, ( 28 ) 
A/t=2ewq (29) 


holds approximately in equilibrium state, by combining Eqs. (26), (27) and (28). Having 
had no report that guesses different peculiar production rate at Kew from those at 
usual land stations, the authors regard the value 4/r as substitutable for Kew. This 


interpretation seems to be supported by the numerical values reported by Chalmers 


(1957) as the characteristic values, that A=3x10-"5 and t<20sec. Thus the assumed 


values may be able to express approximately the practical nature on the value 4/r at 


oop ete calculated results are shown by solid curves in the figure for some field 
intensities at the earth’s surface. 
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Expected Field Profiles 
Solid curves show the profiles in the case that small ions’ life time is taken into consider- 
ation. Dotted curves show the profiles depending on Wilson’s space-charge theory. Curve 
(A) and (a) show the profiles for X,;=30 V/cm, Curve (B) and (b) for X,=50 V/cm, Curve 
(C) and (c) for X,=70 V/cm and Curve (D) and (d) for X,=100 V/cm. 


In comparison with the model proposed here and the results from Wilson’s theory, 
a remarkable difference is found. That is on the increasing manner of field intensity 
with height; the former shows so rapid increase that the intensity at 100m high is 
‘more than ten times as large as that at the earth’s surface, in particular, the increase 
is concentrated within the lowest layer of some ten meters depth, when the surface 
intensity remains less than 30 V/cm, while the latter shows comparatively gradual 
increase, which does not result in the intensity increment by an order of magnitude 
below the levels of several hundred meters high. In fact the profiles for extremely 
intense surface field may lie somewhat toward those of Wilson’s theory with reference 
to concentration of condensation nuclei, some ten thousands per cc at much polluted 


district as Kew. 
3. Reexamination of the Method of Altielectrograph Measurement 


The purpose of reexamination of the Altielectrograph results is generally concerned 
with the relation between point-discharge current in an isolated system and ambient 
field. In addition, the resolution of the trace of the apparatus must be considered 
under the condition in which the rapid increase of field intensity is expected near the 
earth’s surface. The purpose of this section is to point out the importance of the 


resolution. 
The Altielectrograph measurements were essentially based on the similar relation- 


ship to that of Whipple and Scrase, Le. 
i=aX’, ; (30) 


where a=3X10-"2 deduced from the measurement. Accepting Eq. (30), Simpson and 
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Scrase (1937) reported in their paper that ” Similar comparisons between the gradient 
at the ground and the current in the sounding appratus immediately after it leaves 
the ground were attempted, but the number of suitable occasions were too small to 
comfirm the form of the above relationship. If, however, this form is assumed to hold, 


these comparisons indicate the value of @ to be such that 
49 =0.03 F? | approximately, 


when i) is measured in microamperes and F in V/cm.; M can _ be neglected 
when F is not small.” Consider the ascending apparatus, the words “immediately 
after...” can be replaced by the words “floating at levels near the ground.” The 
levels may be decided by the resolution of the apparatus having a discharging tail- 
wire of about 20m long and a recording chart on which the current was recorded by 
the locus of 0.1~1.0mm width and about 15mm length for the first one thousand 
meter. Since the apparatus could not have recorded the current untill it ascended 
more than 20m high, it may have reached at levels of some ten meters high, for the 
reliable comparison with the field intensity at the earth’s surface. To estimate the 
numerical value of coefficient, a, the comparison must be with the measured current 
and the field intensity at a level at which the apparatus is floating. Thus it is probable 
that the apparatus hardly offers the rapid increase of field intensity, and moreover 
makes the comparison tO be nonsense, though it is unimportant with reference to 
Wilson’s formula. It is deduced from the record shown by Simpson and Scrase that 
the field intensity at the earth’s surface remains within 30 V/cm in most cases in which 
the Altielectrograph measurements were held. Then Eq. (31) must amend to 
in=3x10-0(F) (31’) 
1.é. a=3>x10-" 


Eq. (31’) indicates the erroneous estimation of X by an order of magnitude. Thus it 
may be reasonable that the Altielectrograph results showed nearly constant field below 
cloud, against the fact. 


4. Concluding Remarks 


Wilson’s space-charge theory seems to give no practical features because of neg- 
lection of immobilized ions, and moreover, bring the erroneous estimation of the absolute 
field intensity in the Altielectrograph measurement by an order of magnitude. This 
may be able to interpret the gaps between the results from the Altielectrograph and 
Gunn’s observation. The idea that the process immobilizing small ions on large ions 
has an important role for the field profile was also suggested by Adkins (1959). And 
the results depending on it are consistent with the result of Chapman (1958), reporting 
that the numerical value of @ must be corrected to 2.81074, 

To confirm the practical features, it is necessary to measure directly the field 
intensity or ion number, taking careful attention to the rapid increase of field intensity. 
On this point, the observational method similar to the Altielectrograph may be unsuitable 


Electric Field below Cloud 147 


for resolving the discrepancy, but the method that enable to measure directly the ab- 
solute value independent of the field at the earth’s surface is desired. On the ion 
capturing process by rain drop, it may be better to reexamine on the quantitative 
treatment. 
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Abstract ° 


Magnetization-blocking process in chemical remanent magnetization was 
demonstrated through the size development of ferromagnetic precipitates generated 
by an aging of Cu-Co alloy, the magnetic behaviour of which changes from super- 
paramagnetic to stable single domain and then to multidomain one, according to 
the increase in the size of the precipitates. It has been shown that a remanent 
magnetization is acquired during the transition from a superparamagnetic to stable 
ferromagnetic state and that the remanent magnetization thus obtained decreases 
as a result of the formation of domain walls if the size of the precipitates exceeds 
the critical single domain grain size. The result can be well explained by Néel’s 
thermal fluctuation theory of single domain grains and the domain structure theory 
of ferromagnetics. The particle sizes at various stages of aging were estimated by 
different granulometry methods with roughly consistent results. The general idea 
confirmed in the present study has been successfully extended to the CRM phenomena 
in synthetic and natural ferromagnetic minerals. It was suggested that specimens 
produced by low temperature chemical changes can be classified into four groups, 
i.e. the superparamagnetic stage, magnetically-stable single domain stage, the tran- 
sient stage from single domain to multidomain structure and the perfectly developed 
multidomain stage from the view-point of the nucleation of ferromagnetic minerals. 
The magnetic behaviour of several specimens already reported can be consistently 
interpreted on the basis of the present experiment. 


1. Introduction 


It has been shown by Haigh (1958) and the present author (Nagata and Kobayashi 
1958, Kobayashi 1959) that Chemical Remanent Magnetization, or CRM is acquired 
by magnetic materials undergoing a chemical change (e.g. a-Fe,0;>Fe,O0,) in a weak 
magnetic field at a constant temperature below the Curie temperature. The author has 
also shown that the stability of this type of remanent magnetization is very similar 
to that of thermoremanent magnetization (TRM), although the intensity of the former 
is not so great as that of the latter. As for natural rocks and ore deposits, on the 
other hand, many specimens have been presumed as examples of chemical remanent 


magnetization. The intensity of their natural remanent magnetizations is widely 
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distributed. For example, natural specimens containing maghemite or titanomaghemite 
can be classified into two groups with respect to the magnitude of NRM and Q,-ratio; 
the one has a large magnitude of remanence and Q, amounting to 100, while the Q,, of 
the other is quite small and its NRM is sometimes unstable. 

It was first suggested by Haigh that a single domain grain theory of TRM pro- 
posed by Néel (1949) is applicable to the magnetization-blocking process of chemical 
remanent magnetization. The similarity of TRM and CRM in their stability seems to 
indicate the adequacy of Haigh’s theory, although it has not alluded to the variety in 
intensity of CRM at all. 

Experiments demonstrating the production of CRM have recently been carried out 
by the author with Cu-Co alloy (Kobayashi 1960). The alloy has been investigated 
by many researchers (e.g. Mitui and Miyahara 1953, 1955, Mitui 1958, 1959, Sato and 
Mitui 1959, 1960, Weil and Gruner 1956, J.J. Becker 1957, 1958 (a), (b), Bean and Jacobs 
1955), since the phase diagram was established by Tamman and Oelsen (1930) and 
Hasimoto (1937). In the alloy, cobalt is soluble only to the extent of 5% at 1100°C, 
3% at 940°C, and 0.3% at 400°C. “It is, therefore, possible to obtain the fine ferro- 
magnetic precipitates of face-centred cubic structure containing approximately 10 % 
Cu in the non-ferromagnetic Cu matrix, if the homogeneous Cu-Co alloy containing 
less than 5 % Co is annealed at a temperature of about 600~800°C. 

Since the Curie temperature of the precipitates from the Cu-Co alloy is far above 

the corresponding solvus curve, variation of magnetization and remanent magnetization 
of the precipitates with increasing aging time can directly been measured at the aging 
_ temperature. It has been concluded from the varied phenomena observed by other 
authors that the precipitation proceeds in this alloy in such a way that the average 
particle size of the ferromagnetic precipitates increases with time, while the total 
amount of the precipitates is kept constant during aging except for a short time interval 
in the initial stage. From these view-points, this alloy may be the most suitable specimen 
to examine the magnetization-blocking process during chemical changes. 

Measurements carried out so far by the author are as follows: 

(i) Change in saturation magnetization with aging time. 

(ii) Total and remanent magnetization at the various stages of aging under the 
influence of various intensities of a magnetic field. 

(iii) Production of thermoremanent magnetization by cooling the precipitates under 
the influence of a magnetic field. 

The phenomena observed in (ii) must be distinguished from the so called magnetic 
field annealing effect. The former is related only to the alignment of the magnetic moment 
of each particle, while the latter concerns change in the bulk properties such as ani- 
sotropy, coercive force and so on, induced in the specimen. 

In the following sections, it will first be mentioned how to prepare and how to 
heat-treat the specimen (§ 2), then the results of measurements (i) and (ii) will be described 

in § 3. In § 4, a theoretical interpretation of the measurement (ii) will be proposed. 
Detailed and quantitative discussion of the results is described in § 5. § 6 treats (iii), 
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in which the size distribution of superparamagnetic precipitates is inferred from the 
partial TRM characteristic of the specimens. In § 7, the theory derived from the 
present experiments on Cu-Co alloy will be consistently extended to the interpretation 


of CRM in both natural and synthetic specimens. 


2. Preparation of the Specimens and Heat Treatment 


Alloy specimens were made from electrolytic copper of 99.9 % purity and cobalt 
of 99.8% purity. A little more than 2 weight per cent cobalt and balance copper 
were mixed and melted in an induction furnace in a vacuum of nearly 10-*mm Hg. 
The ingot of the alloy was cut into rods of about 4cm in length and nearly 0.8cm 
in diameter. Sealed into a quartz tube, it was homogenized by heating at 1000°C for 
more than 50 hours in a hydrogen atmosphere. After the heat treatment, it was rapidly 
quenched with water by crushing out the sealing tube. The specimens thus prepared 
are paramagnetic or weakly superparamagnetic, their susceptibility being less than 10~* 
emu/gr. Heat treatment for the formation of a uniform solid solution from the unmixed 
specimens, after an aging experiment was usually completed within several hours at 
ca. 1000°C. 

In the case of the aging experiment, the specimens of uniform Cu-Co alloy, sealed in 
an evacuated quartz tube, were inserted into the non-inductive electric furnace after 
the temperature in the furnace had settled at the designated temperature (usually 750°C). 
The specimens were expected to arrive at the designated value of temperature within 
a few minutes. Then the specimens were aged at the constant temperature. 


3. Changes in Magnetic Behaviour with Aging Time 


3. 1. Measurement of saturation magnetization 

Fig. 1 shows the magnetization curves towards saturation of the five kinds of 
specimens after quenched in various stages of aging at 750°C to room temperature. 
For all kinds of the aging time intensity of saturation is practically constant. It may, 
therefore, be safely said for the present specimens that precipitation of the Co-rich 
phase is finished within 5 minutes, and afterwards, the change in the observed proper- 
ties is associated exclusively with the increase in average particle size of a constant 
amount of a precipitated ferromagnetic phase. 


3. 2. Total, remanent and reversible magnetization at various stage of aging 

Magnetization of the alloy at various stages of aging under the influence of various 
intensities of a magnetic field was measured by means of a high temperature ballistic 
method designed by Uyeda (Nagata, Uyeda and Akimoto 1952). Remanent magnetization 
was also measured with the same method by switching off the electric current of the 
field coil for several seconds. The influence of the suppression of the applied field for 
such a short interval may be negligible in the results. Fig. 2 (1), (2) shows two examples 
of the results of these measurements. 

Reversible magnetization /,, mostly due to superparamagnetism (except in the later 
stages of precipitation) was calculated here as the difference between magnetization . 
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Fig. 1. Magnetization curves of Cu-Co alloy after quenched from various 
stages of aging (measured at room temperature) 


a02 


000 ‘ —— 0! 
1 2 3 4567870 20 40 60 400 200 400 600 1000min 07 1 10 00 min 
mie sates —+ t\log -scale) 
(Maoh Hel ‘Os, Q\ioih iHep=20 Oc. 


—Q— total magnetization J 

—@— remanent magnetization /, 

-~--- reversible magnetization Jp 

Fig. 2. Total, remanent and reversible magnetization of Cu-Co alloy 
at various stages of aging at 750°C. (measured at 750°C) 


with the applied field J and remanent magnetization /,, namely, from the eqation: 


J=J-+Jo- 


It must be noted in the figure that the intensities of reversible magnetization and 
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remanent magnetization arrive at their maximum values after the constant time interval 
of aging, say nearly 10 min and 50 min respectively, almost independently of the 


intensity of the applied field. 


4. General Aspects of the Phenomena 


4. 1. Superparamagnetic stage of precipitation 

It is well-known (Kittel 1946, 1949, Néel 1947) that a ferromagnetic particle changes 
its magnetic structure from multidomain to single domain when its size becomes quite 
small. Further, Néel (1949) pointed out that if a single domain particle is small enough, 
its direction of magnetization may be spontaneously reversed by the effect of thermal 
activation energy k7' which overwhelms the energy barrier against the rotation of 
magnetization K,v even in the absence of a magnetic field. 

The process is characterized by the relaxation time tp), that is e say, the time 
decay of a remanent magnetization of an assembly of such particles is expressed by 


JAt)=J0) exp (—t/t9) (1) 
The relaxation time ty is given as 
1/to=fo exp (—Kywv/kT), (2) 


In this expression, v is the volume of the particle, f is a frequency factor of the 
order of 10° sec™!, and it is assumed that the particles have a uniaxial anisotropy the 
first order coefficient of which is K,v. For particles with an anisotropy of cubic sym- 
metry, the energy barrier between the adjacient easy directions, ie. —Kv/4kT for K>0 
and —Kv/12kT for K<0, will appear in the exponential instead of K,v in the above 
formula (2). 

The above treatment is for the case of a zero applied field. If a reversing field H 
is applied along the easy axis of an aligned assembly of uniaxial particles, the energy 
barrier for rotation of magnetization is K,v (1+/,H/2K,)*, or approximately K,v—vjJ,H 
for a small H. Then the equation (2) is changed as denoted by 


V/ty=fol(1 + J,H/2Ky)(1—J2H?/4K,2)'exp (—(Kyv/RT\(1+J.H/2K,)°) +(1—J,H/2K,) 
(1—J,2 H?/4K,2)*exp {—(K,v/kT)(1—J,H/2K,)"} ] (as 


Because of the exponential dependence of r (rt) or t,,) on v/T, there is a fairly 
well defined value of a,=(v/T), at which the transition from t<i, to r>#, occurs, 
where ?¢, is the time for experiment. When v/T is less than a,, the assembly can come 
to a thermal equilibrium state in the experimental time. This state can be achieved 
either when the volume of each particle v becomes smaller than v,=a,T at the ordinary 
temperature or when the temperature is higher than T,=v/a,, which is called “blocking 
temperature”. Table I shows the calculated particle radius for various materials corre- 
sponding to t)=10* sec at room temperature, assuming either one of the two kinds of 


origins of anisotropy, ie. magnetocrystalline anisotropy K (shape is spherical) and 
maximum shape anisotropy z/,.. 


——— 
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Table I. Critical radius corresponding to ty=10?sec for various materials (at room temperature) 


Origin of anisotropy Fe spe wa Fe,O, | 7-Fe.0; | a-Fe,0, 


Magnetocrystalline A 
anisotropy K, or Ky; 130 40* 140* 290 385 ~ 290A 


Shape anisotropy** 
aJ,? 


30* 35 35 70 85 2240 A 


* After Bean and Livingston (1959) 
** In this column the radius of the equivalent sphere having a volume equal to the critical 


one is listed. 

In a thermal equilibrium state, particles behave magnetically in a manner analogous 
to the Langevin paramagnetism of moment-bearing single atoms. As has been well 
known, the magnetization curve of the assembly of identical particles with volume v 
can be denoted by the familiar Langevin function with quantum number j=o, ic. 


J/J,=L.(a)=cotha—1/a (4) 
where 


a=vj],H/kT 


provided that the anisotropy of the particles is assumed to be negligible compared 
with the thermal agitation. For a low field this function is approximately given as 


J/Js~a/3=v],H/3kT (5) 
For a high field, it approaches saturation in the form of 
I/Je~A-1/a=1-(kT/0J.)* Fy (6) 


This behaviour of small particles has been termed “ superparamagnetism ” by Bean 
(1955). 

In the present aging experiment, ferromagnetic precipitates are so small in particle 
size that they behave “superparamagnetically” in the first stage of aging. Thus, 
according to the increase in the average particle volume v with increasing time, re- 
versible magnetization of the system in a low magnetic field increases with aging time 
as clearly seen in eq. (5) in which ///, is nearly proportional to »v. 


4. 2. Stable single domain stage 

When v/T exceeds a, according to the development of the volume of the precipi- 
tates, relaxation time is much larger than the experimental time and the magnetization 
of particles becomes stable. The state of magnetization of the system is then “ frozen” 
and the specimen will maintain a memory of the field present at the time of v/T=«,. 
Fig. 2 indicates these features quite well, in which remanent magnetization developes 
abruptly within short time intervals after the appearence of the maximum super- 
paramagnetism. 

The freezing of magnetization may occur on account of at least two kinds of 
- mechanisms; one is the cooling of the system from a temperature above its blocking 
temperature 7;=v/a,, which was treated by Néel (1949, 1955) in his theory of the thermo- 
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remanent magnetization of single domain particles, and the other is the case of the 


generation of remanent magnetization treated in the present study. 


4. 3. Multidomain particle stage 
Another interesting feature of Fig. 2 is the gradual decrease of remanent magnet- 


ization in the later stage of precipitation. It may be interpreted as a result of the formation 
of magnetic domain walls in precipitated particles. In a sufficiently developed particle, 
it may be easily understood that memory of the direction and intensity of magnetization 
obtained in the single domain stage can scarcely been maintained because of the for- 
mation of a large number of domains, directions of which are distributed in such a way 
as to minimize the total magnetic energy. The transition stage from single domain to 


TH 
0-O-U- 


(Q) (b) (Cc) (d) (e) 


(a). superparamagnetic stage 

(b). magnetically-stable single domain stage 

(c). transition stage from single domain to multidomain structure 
(d). multidomain without surface closure domains 

(e). multidomain stage 

Fig. 3. Schematic illustration of the process of the formation of 
the reversing domains in uniaxial crystal with increasing volume. 
This domain structure must be slightly modified for cubic crystal. 


multidomain has scarcely been studied as yet. Fig. 3 shows a schematic view of this 
transition of domain structures of uniaxial material such as MnBi. The situation seems 
to be similar in Co precipitates with cubic anisotropy. 


5. Detailed Discussions 


5. 1. Field dependence of peak intensity of the superparamagnetic component. 

As already mentioned in section 3, curves of the change in intensity of the super- 
paramagnetic component in magnetization with aging time have a peak nearly at a definite 
position in t-axis, say nearly 10 min, independently of the applied magnetic field. This 
result indicates that the peak appears at the time when a definite average volume vy 
of precipitates is realized, since the mode of precipitation is the same in every experi- 
ment. So that it is possible to discuss the field dependence of the magnetization of a 
superparamagnetic assembly consisting of particles with the average volume vo, by 
plotting the peak intensity of such curves as shown in Fig. 2 as a function of applied 
field. Fig. 4 shows the field dependence of the superparamagnetism (hollow circles) 
thus plotted up to 500 Oe. 

Considering an assembly consisting of identical particles with volume vo, the /~H 
relationship of the system can be denoted by the Langevin function L..(@) as already 
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Fig. 4. Field dependence of peak intensities of superparamagnetism 

and remanent magnetization 
introduced in 4.1. In the actual case, however, the relationship must be modified a little 
because the volume of particles is not exactly equal to v, but is widely distributed 
around vy. Kneller (Kneller 1958, Hahn and Kneller 1958) successfully treated this by 
assuming a rectangular distribution of the volume. 
Now let us assume, according to him, that the volume v is uniformly distributed 
between v,(1—d) and v,(1+4). Then the formula for magnetization is given as 


(+b) (1—b)sinh[ao(1+0)]) _ 
w= (1/2bay)fn" , Lala)da=(1/2ba,) In| of = 1 7 
Jo/Tn=(h/2bay)\ i y Pala da=(1/2bae) Inyo py enae) (7) 
where a,=v)/;H/kT and /,, /p,, denote respectively the intensity of the saturation magnet- 
ization of a precipitated phase and the saturated value of the superparamagnetism of 
the alloy. The above function L*(a) is called a “ modified Langevin function ”. 


Since D=v, and (1/v)=(1/250,) In +? (8) 
an approximated representation for low fields (a<¢1) is 
T/T e~ 0/3 = V0JH/3kT (9) 
and for high fields 
JolJou~1—CJa)=1—(1/2ba,)|nf ++ | (10) 


In other words, this sort of distribution of the particle volume affects the magnetization 


curve, only near the approach to saturation. 

The upper curve in Fig. 4 is a modified Langevin function thus obtained with 
Vp=3.8 X 10718 cm’, (v>=100A), b=0.71 and Jy,,=14emu/cc. In this calculation, the 
spontaneous magnetization of a precipitated substance (~90% Co-Cu) is taken as 
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J,=1300 emu/cc. The curve is fairly well fitted to the measured values (hollow circles). 
This result indicates that the precipitated particles in the specimen examined here are 


between 66A to 120A in their radii. 


5. 2. Field dependence of the peak intensity of remanent magnetization. 

Although the curves of remanent magnetization as a function of aging time also 
reach a peak after a constant aging time, a mathematical treatment of the field dependence 
of their peak values is more complicated than in the case of superparamagnetism. 
For the sake of simplicity, consider an ideal model of the alloy in which N number 
of identical particles with equal volume v precipitated in non-magnetic matrix. With 
the increase in aging time, the volume of each particle v increases, while N decreases 
in such a way as to keep the total volume of the precipitated materials Nu=V always 
constant through the precipitation process. 

Fig. 5 shows schematic curves representing changes in the reversible and remanent 
magnetization with particle volume v at a constant temperature under the influence of 
a magnetic field H far from the saturation field. In the abscissa of the figure, v, and 
v,’ denote the critical volumes corresponding to v/T=a, in a zero field (from (2)) and in 
H=H, (from (3)) respectively. From eq. (3), ve’ can be approximated by vu o(1— J,H/2Ku)-? 
for uniaxial particles. If the symmetry is not uniaxial, it is sufficient to replace Ky 
by CK in the above parenthesis, where C is a numerical constant dependent on the 
symmetry. 


Fig. 5. A schematic representation of change in re- 
versible magnetization (upper curve) and remanent 
magnetization (lower curve) with particle volume in 
a simple model. ve~vp, J4~Jg are given in the text. 


Assuming that the particles have a spherical shape and an infinite number of 
easy magnetization axes, the intensities of J, and J, are given as follows: 
1. v<v,; superparamagnetic stage 


Jp=Ja=(V/U J Le(vJ,H/RT) (11) 
since Jp,,=(V/U)/,, where U is the valume of the alloy specimen. 
2. v=v,; the beginning of magnetization-freezing, 


J,-=Jp=(V/U)J Ll VeJsHo/kT) , (12 ) 
Se ee Says eg < ke zE 


a ow 
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Jr=Jo=(V/U VJ LevJsHo/kT),  Iv=Sn (13) 
4, vVp=v>v,'; the magnetic moments are completely frozen, 


J-=Jn=(V/U Le(veJeHo/kT), Jn=Sn (14) 
5. U>Vp; formation of domain walls 


Jo>Jr>0, Jn>STe 


where J;/H) shows the ordinary reversible susceptibility of stable ferromagnetics. It 
must be noted that, as seen in the above equation (14), the peak intensity of remanent 
magnetization is dependent on »v,’, not on vp. 

In the actual alloy, the precipitates have various sizes with different rates of size 
development, and moreover, they have f.c.c. crystal symmetry and, in a later stage of 
precipitation, a plate-like shape (Bean, Livingston and Rodbell 1957). Nevertheless, these 
formulae seem to be applicable to the analysis of the measured results. 

Néel assumed, in his theory of TRM, the assembly of prolate-spheroid particles, 
1/3 of which have major axes parallel to the field and 2/3 of which have those 
perpendicular to the field. According to his model, the above mentioned formulae (12), 
(13), (14) must be replaced by 


J,=Jz=(V/U \(J./3)tanh(vJ,Ao/kT) (15) 


etc. In the present case, however, L..(a) seems to be preferable because the precipitates 
in the alloy are oriented at random without any preferential directions. 
The lower curve in Fig. 4 illustrates the Langevin function L..(v¢J;H)/kT) estimated 
from the numerical value of the critical volume given so as to fit the curve to the 
measured values (solid circles) as closely as possible. It-is remarkable that the curve ap- 
proaches saturation in an applied field below 100 Oe, which is very weak compared 
with the coercive force of the same alloy already aged to a magnetically-stable stage. 
From this curve v, is determined to be 8.8x107!%cm’, then y,=128A. This value of 
r. seems to be in agreement with one estimated from v/T=a,, T=750°C, which is 
about 130A if an approvable spheroidal shape with the axis ratio 1.4 (M—N=0.22) is 
assumed for the precipitates. This is also quite acceptable since it is a little larger 
than r determined from the curve of maximum superparamagnetism. 


6. Thermoremanent Magnetization and Size Distribution of Superparamagnetic 
Particles in Various Stages of Aging 


As mentioned in the previous section (3. 2), magnetization of superparamagnetic 
particles is blocked by cooling through their blocking temperatures as well as by the 
increase in the particle size. Fig. 6 shows the growth of remanent magnetization in 
the case of cooling appropriately annealed Cu-Co alloys to a low temperature in a 
magnetic field of 10 Oe. Measurements are performed as follows: at a designated 
stage of aging at 750°C the specimen is cooled down to the liquid MN, temperature 
under the influence of a magnetic field (10 Oe), the thermal variation of remanent 
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Fig. 6. Growth of TRM in Cu-Co alloy after various intervals 
of aging at 750°C. 
magnetization being measured at various temperatures by means of the ballistic method. 
Since the precipitation does not proceed at temperatures below approximately 500°C, 
the remanent magnetization thus acquired in the specimen is a typical example of 
“thermoremanent magnetization” (TRM) obtained by each precipitated particle when 
it is cooled through its blocking temperature. It is also shown by experiment that, 
if the specimen is heated, the remanent magnetization inversely decreases with in- 
creasing temperature along the same curve as in the cooling process. This is a well- 
known property of TRM. 
In an assembly of particles having different blocking temperatures of a wide range 


such as the precipitates in Cu-Co alloy, there holds the linear addition law of partial 
TRM, namely, 


n Ti-y To 
pH Tin (To)=J 7 (70) ( 16 ) 


where / a (1) denotes partial TRM, i.e. a remanent magnetization after cooling to 
To, H being applied only in a temperature range between 7, and deen (I> a 
the saturation magnetization of f.c.c. Co-precipitates is practically constant between 
T, and the aging temperature T,,, J tte may be replaced by / phaelie in the present 
case, and thus partial TRM for 47=7;—T,_,=100°C can be estimated as 

1b ml ee at various values of 7; from the curves in Fig. 6. 


Fig. 7 demonstrates “spectra” of the partial TRM thus estimated. Since blocking 
of magnetization occurs at v/T=o, remanent magnetization /J,/J,,,= ut ai ae is 
independent of the particle volume. Therefore, it is concluded that the height of each 
column of partial TRM in the spectra is proportional to the relative content of particles 


having blocking temperatures within the corresponding temperature range. Particle 
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Fig. 7. Spectra of partial TRM produced in different 

temperature ranges in appropriately annealed Cu-Co 

alloy. 
radius inferred from v/7;=a, is also shown in the abscissa of the figure, assuming 
r=128A for 7,=750°C. The average radius of precipitates in the specimen aged for 
12 min is estimated to be 117A. This value is a little larger than 7=100A which is 
calculated from the field dependence curve of the maximum superparamagnetism. The 
possible reason for this discrepancy is that the precipitation proceeds for nearly 5 min 
further before the temperature of the specimen falls below 500°C. 

On the right-hand side of the figure, the relative intensity of remanent magneti- 
zation blocked by volume development during aging which is read from the heights 
of the right ends of the curves in Fig. 6 is also demonstrated. It can be clearly seen 
from the mode of partial TRM spectra how the distribution of the size of the pre- 
cipitates gradually shifts towards the right side with increasing aging time until the 
magnetization of particles is frozen one after another from larger ones during an 
isothermal aging. 

It must be mentioned that such spectrum of partial TRM as those described here 
may be useful as a technique for determining the distribution function of particle size. 
If the saturation magnetization of the fine particles /, is not constant in the whole 
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: T : : 
range of temperatures, it is necessary to use the ratio 4/ Pipe ‘)/Js(T;) in place of 


aye 


Tn : 
method. Methods of determining the distribution of particle volume have been proposed 


by several authors (Weil 1954, Weil and Gruner 1956, Henning and Vogt 1957). Weil 
proposed a method to estimate the size distribution of superparamagnetic particles 


(T,) as the measure of the relative content of particles in the granulometry 
»H 


from the temperature dependence of the ratio of remanence to saturation magnetization. 
The method proposed in this section is similar to Weil’s but seems to be more readily 
available to us. 


7. Interpretation of Chemical Remanent Magnetization on the 
Basis of the Cu-Co Experiment 


It seems that the comparatively simple treatment of the experimental results per- 
formed in the previous sections is based on the facts that the precipitation in Cu-Co 
alloy is the simplest process of the volume development of precipitates, and that the 
magnetic interaction among the particles is negligibly small because of the dilute 
distribution of the precipitates in the alloy. In the case of general chemical reactions 
such as reduction of a-Fe,0; to Fe,O, and oxidation of Fe,0, to y-Fe,O,, etc., the 
problem is more complicated and may be difficult to examine with exactness in experi- 
mental processes. Nevertheless, the similarity in the stability of CRM and TRM seems 
to support the possibility of extending the conclusion obtained in the case of Cu-Co alloy 
for interpreting CRM in general cases. It will be shown in the following that magnetic 
properties of both natural and synthetic specimens can be successfully explained on 
the basis of the Cu-Co experiment. 

It may be possible to consider a chemical change as a kind of nucleation process. 
Considering the case of the reduction of a-Fe,O; to Fe,O,, for example, a very small 
amount of Fe;O, will be formed probably around a dislocation or other imperfection in 
the a-Fe,O, grain in the early stage of the reduction. Such nuclei are superparamagnetic 
and can be easily magnetized with the applied field on account of their fine sizes. As 
the chemical reaction proceeds, these nuclei will grow, and rapidly become magnetically- 
stable and their magnetization will be frozen along the direction of the ambient field. 
After that, the reaction will still continue to proceed until the whole volume of the 
grain is completely filled with Fe,Q,. 

If the final size of the nucleus is smaller than the critical size for a single domain 
grain, blocked magnetization will be observed as a CRM. On the other hand, if the 
size is larger than the critical one, the nucleus becoming multidomain in structure, 
magnetization once frozen will decrease owing to the formation of the reversely mag- 
netized domains in the same nucleus in a way similar to the Cu-Co alloy. In the 
extreme case where the nucleus grows to a considerably large size compared with the 
critical size, CRM will disappear so that the specimen will have only the isothermal 
remanent magnetization in the applied field. 


In the experiment of the generation of CRM by chemical changes such as the reduction 
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of a-Fe,O, to Fe;0, and the oxidation of Fe,;O, to y-Fe,O; already reported by the 
present author (ibid.), it was confirmed by an electron-microscope observation that the 
grains of a-Fe,O; and Fe,O, have a roughly spherical shape with diameters of 0.1 u~ 
0.6. Since the critical single domain size of Fe,O, and 7-Fe,O, with spherical shape 
is about 0.14:~0.2 according to Morrish and Yu’s estimation (Morrish and Yu 1955), 
_ most of the grains in these powders are multidomain in structure. Thus, it can be 
said that the CRM of these specimens should be reduced by the domain wall formation 
compared with the case when the grain size is a single domain one. This conclusion is 
in general agreement with the experimental result indicating that the intensity of the 
CRM of Fe;O, is smaller than that of the TRM by nearly one order of magnitude when 
the applied field is small. 

In the case of oxidation of Fe,O, to y-Fe,O;, it has been ascertained that the intensity 
of CRM is nearly twice that generated by the reduction of a-Fe,O; to Fe,;O,. This 
fact seems to be quite reasonable because the estimated critical single domain grain 
size of y-Fe,O; is nearly twice that of Fe,Q,. 


8. Various Phenomena in Natural and Synthetic Specimens 
releyant to Superparamagnetism and CRM 


From the stand-point of the same theoretical idea as mentioned in the present 
article, the magnetic behaviour of various kinds of natural and synthetic specimens can 
be consistently explained. Characteristics of several specimens already published are 
listed in Table II. The stages of nucleation are classified here roughly into four, namely, 

_the superparamagnetic stage, magnetically-stable single domain stage, the transient 
stage from single domain to multidomain structure and the perfectly developed multi- 
domain stage. et 

Each specimen can be successfully classified into each stage of nucleation according 


to the particle size of the concerned mineral. 


9. Concluding Remarks 


In this article the magnetization-blocking process in chemical remanent magneti- 
zation and the variation in the intensity of remanent magnetization due to domain 
wall formation according to the size development of the precipitated ferromagnetic 
materials were clearly demonstrated with Cu-Co alloy. It was confirmed that the result, 
being fairly well explained by Néel’s thermal fluctuation theory of single domain 
grain and the domain theory of ferromagnetics, can be successfully extended to the 
various kinds of CRM phenomena. 

However, some of the important problems in CRM still remain unsolved. Among 
these, the following two seem to be the most important, i.e. 

(1). the magnetic interaction between the precipitates and matrix substances as 
well as among the precipitates themselves. 

(2). the effect of the magnetic anisotropy of ferromagnetic precipitates on CRM. 


Studies on these problems are now being undertaken. 
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_ Abstracts of the Papers Presented at the Symposium on Rock Magnetism, 
at Shinshu University, Matsumoto, Oct. 5-6, 1960. 


SYONO Y. 
University) 


(Geophysical Institute, Tokyo 
Remanent Magnetization of 
Large Single Crystals of Ferromagnetic 
Minerals Field dependence of TRM and 
IRM and their stability against AC demagneti- 
zation of ferromagnetic single crystals of natural 
hematite, phrrhotite, magnetite as well as 
synthesized YFeO;, magnesioferrite, yttrium 
iron garnet were studied. It was confirmed 
that characteristic features of TRM still remain 
in fairly large single crystals. The development 
of TRM is effective in the crystals having 
strong magnetic anisotropy, especially in para- 
sitic ferromagnetic crystals such as YFeO, and 
hema tite, while in the cubic crystals the value 
of Jy-/J/s which may represent a degree of fixing 
of TRM is only 107° at most but even in this 
case TRM can bé definitely distinguished from 
IRM. TRM of the former is far more stable 
than that of the latter against AC demagneti- 
zation. These experimental results may be 
explained as follows, the larger the magneto- 
crystalline anisotropy is and the smaller the 
spontaneous magnetization is, the larger the 
coercive force and the more remarkable the 
- fixing of TRM becomes. 


AKIMOTO S. and YAMA-AI M. _ (Geo- 
physical Institute, Tokyo University) Imper- 


fect Reverse TRM of II-Ht Series 
TRM of II-Ht series (Il: 0~50%) was examined 
in detail with special reference to the production 
of imperfect reverse TRM which was found 
by Uyeda and Carmichael in a certain compo- 
sitional range of the II-Ht series (I! 10~15%). 
Ten specimens were newly prepared by Akimoto’s 
method, their composition determined by X-ray 
analysis being (5), 7, 14, 20, 25, 27, 31, 35, 40, 
449 Il. All ten specimens show more or less 
a character of imperfect reverse TRM in the 
original state (1150°C in vacuum, 6 hrs, quench); 
for example, 44% Il specimen’s TRM decreases 
with decreasing temperature in the temperature 
range from room temperature to —196°C and 
a reversal temperature at which the self-re- 
versal of TRM should take place was extra- 
polated to be —240°C. The effect of heat and 
pressure treatment on this reversal temperature 
was examined. Reversal temperature of the 
specimens with 44~27% was remarkably ele- 


vated by annealing at 500°C for 112 hrs.; TRM 
of the specimens with 40 and 44 I1% was 
reversed at room temperature. Three specimens 
with 35, 40 and 44 I1% were also subjected 
to pressure of 2500 bar for 50 hrs. at 300°C 
but the effect of such pressure treatment was 
found to be far small compared with that of 
heat treatment. 


(Geophysical Institute, 
Chemical Remanent 
Magnetization VII (Superparamagnetism 
of 7 Fe,0; Powder)——Superparamagnetic 
fine powder of 7 Fe,O; was artificially produced 
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_ by dehydration of yFeOOH. As X-ray analysis 


shows the co-existence of aFe,O;, 7Fe,O; is 
finely dispersed in particles of aFe,O,. Mag- 
netization curve of the specimen is completely 
reversible with H, and approaches saturation 
even in the applied field of several hundreds 
Oe. It can be quite well represented by modi- 
fied Langevin function if the particle radius of 
7 Fe,O,; is assumed to be between 80 A and 290A. 
This specimen is transformed to stable ferro- 
magnetic materials when cooled to about 80°K 
(liq. N. temperature). If it is cooled in a mag- 
netic field, TRM is produced. The blocking 
temperature of the 7 Fe,0, powder is determined 
to be —140°C at most by measuring the thermal 
decay of TRM, the result being well consistent 
with that inferred from particle radius above 
estimated. The relation between superpara- 
magnetism and CRM were also discussed from 
theoretical view-point. 


Suimizu Y. (Geophysical Institute, Tokyo 
University) Magnetic Viscosity in Rema- 
nent Magnetization of Magnetites in 
Rocks Time change in remanent magneti- 
zation of magnetite separated from iron sand 
and some sedimentary rocks was measured 
after the suppression of magnetic field. As 
already reported, remanent magnetization Jr 
varies with time as 

Jr=Jro—S (logt—log ?t’) 
Thus, the magnitude of S was determined for 
various specimens as a function of their coercive 
force and magnetic field applied to produce the 
remanent magnetization. The results are fairly 


(165 ) 
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consistent with those obtained in the previous 
studies. 


DOMEN H. (Physical Institute, Faculty of 
Education, Yamaguchi University) Effects of 


the Axial Pressure on Remanent Mag- 
netization For the purpose of under- 
standing the effects of the axial pressure on 
remanent magnetization of rocks, time change 
in remanent magnetization under the influence 
of a small constant axial pressure has been 
examined by using an artificial rock made of 
white cement in which fine particles of magnetite 
are dispersed as uniformly as possible. The 
specimen is formed into a cylinder, 14mm in 
diameter and 10mm in height. It contains 
10 wt% magnetite powder (iron sand of 99% 
purity) collected from Kihado Tertiary sedi- 
mentary bed at Kitaura, Yamaguchi Pref., size 
of which is about 0.3mm and the Curie tempera- 
ture of which is above 280°C. Experiments 
were carried out as follows for both IRM 
acquired at the time of solidification of cement 
in geomagnetic field and TRM produced by 
geomagnetic-field-cooling from above 600°C; 
the specimen is mounted on the press in such 
a way that the vector of remanent magnetization 
is nearly perpendicular to the direction of the 
earth’s field and the axial pressure of 0, 2, 7, 
12, or 30kg/cm? is applied along the vertical 
axis of the specimen for nearly 2,500hrs in 
total. Time change in remanent magnetization 
is measured by an astatic magnetometer in the 
pressure-free state after each designated time 
interval. It was shown by this experiment 
that the horizontal component of both IRM 
and TRM increases with increasing time along 
the external (geomagnetic) field, but that the 
tendency decreases with increase of the pressure. 
These phenomena are more remarkable for 
IRM than for TRM, indicating that IRM is 
more unstable than TRM. On the other hand, 
the increase of remanent magnetization im- 
mediately after the application of the pressure 
(probably owing to PRM) seems to become 
large as the pressure increases. 


SASAJIMA §. (Geological Laboratory, Fukui 
University) Rock-magnetic Considerations 
of Granitic Rocks Magnetic character- 
istics of metamorphic rocks and granitic rocks, 
younger and older, occurring in various meta- 
morphic regions of Japan have been examined 
with special reference to the Hida metamorphic 
complex. The principal results of the study 
are as follows: (1) Metamorphic rocks show 


generally no detectable ferromganetic behaviour 
on thermomagnetic analysis and their NRM is 
less than 10-‘emu/gr. (2) The older intrusive 
granites show very feeble ferromagnetism as 
comparable to paramagnetism, their NRM being 
much scattered in both direction and intensity 
even in the same rock-mass. (3) The younger 
granites have much more ferromagnetic miner- 
als, their NRM being fairly concentrated. (4) 
The petrochemical consideration of the granites 
related to FeO-Fe,O,-TiO, ternary system indi- 
cates that the ratio O/Fe becomes small according 
to the development of metamorphism and that 
the younger and older granites in the same 
metamorphic region are much connected with 
each other in their composition, the former 
being always more oxidized than the latter. 
(5) The above-mentioned results can be explained 
if it is assumed that the position of Fe- and 
Ti-atoms in the rocks is uniquely determined 
by the oxygen partial pressure (i.e. H,O-vapour 
pressure) under which those rocks are generated. 
Namely, if the total rock pressure is constant, 
Fe-atoms take their seats in silicate minerals 
as ferrous ions under very low O,-pressure 
(metamorphic rocks), while under high O,- 
pressure they crystallize as ferromagnetic oxide 
minerals the composition of which is determined 
by the O,-pressure (the younger granites). This 
interpretation will be further developed by 
taking Kawai’s P-T stability relation into ac- 
count. 


YUKUTAKE T. (Earthquake Research 
Institute, Tokyo University) Palaeomagnetic 
Study on the Volcanic Rocks on Osima 
Island Variation of the geomagnetic field 
for the past one thousand years has been 
inferred from the direction of the remanent 
magnetization of rocks ejected at periods known 
by old documents. The results show an 
agreement with those of instrumental obser- 
vations which have been conducted since the 
sixteenth century. Periodical changes are found 
both in declination and in inclination: Incli- 
nation varies rather slowly with the period of 
about 600-700 years, while declination changes 
with 400-500 year period. 


ASAMI E, (Faculty of Science, Kyoto 


University) NRM of Uyamazaki-Basalt, 
Yamaguchi Pref. The palaeomagnetic 
measurement in Yamaguchi Pref., made in 
1950, indicates that among 10 specimens taken 
from Uyamazaki basalt, which is petrologically 
similar to Cape Kawajiri basalt with mixed 
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reverse and normal NRM, there are 8 normal 
and 2 reversed ones. To acertain this result, 
more than 80 basalt specimens were collected 
at outcrops of about 1.5km wide along the 
coast with the interval of nearly 10m to each 
other. It has been shown from the measure- 
ment of their NRM that all of the 80 specimens 
are normally magnetized along the present 
geomagnetic field contrary to expectation. In- 
tensity of NRM, on the other hand, varies from 
0.8x10-* to 24.8x10-°. This fact seems to 
indicate that the possible occurrence of reverse 
NRM might be found if the collection is per- 
formed more in detail, such as with the inter- 
vals of several meters. 


MoMosE K. (Faculty of Liberal Arts 
and Science, Shinshu University, Matsumoto) 
Intermediate palaeomagnetic directions 
retained in rocks For the purpose of 
obtaining detail knowledges about the stability 
of geomagnetic field during Tertiary times, the 
author have engaged in the measurements of 
palaeomagnetisms of rock samples taken from 
the Komoro, the Shigarami and the Enrei for- 
mations in Central Japan. Among the collected 
samples, some of the samples of which the age 
is referred to the middle Pliocene, indicate their 
pole positions calculated from virtually mag- 
netized directions of samples being located in 
- the low latitudes (30°S-30°N) on both sides of 
the present equator. The fact will signify that 
the geomagnetic field might have changed its 
polarity from the time of normal direction (the 
lower Pliocene) downward to the time of reverse 
direction (the upper Pliocene). Variations of 
intensity of the main field through this process 
of reversal, were inferred from the values of 
Jn/Jre obtained from rock samples. The results 
indicate that decrease of geomagnetic intensity 
is assumed for many samples when the pole 
positions were near to the present equator, 
whereas the intensities seem to have taken in 
general the values approximately the same as 
the present, when the pole positions were in 
the high latitudes. The rock samples here 
treated with, were carefully examined to verify 
that their magnetic properties can be fit for 
palaeomagnetic purpose. Examinations of rock 
magnetism consist of (1) A-C demagnetization 
(2) measurements of Js(T) (3) measurements 
In/Jre (4) X-ray analysis and (5) thermal de- 
magnetization. The results of examination 
proved that the ferromagnetic minerals contained 
in these rocks were identifiable wholly with the 
ensemble of solid solutions belonging to the 


ulvospinel-magnetite series, so that they are 
considered to retain stable NRM. Based upon 
the results thus obtained, it seems likely to be 
said that when dipole positions were in the 
low latitudes, the geomagnetic intensities might 
have decreased down to be within the values 
ranging from one-third to even one-tenth of 
that of the present. Dipole positions could be 
obtained simply by calculation from virtual 
magnetic directions of rock samples even when 
the pole positions are assumed to have located 
in the low latitudes. Seeing that the decrease 
in geomagnetic intensity might have taken place 
through the process of reversal, the author is 
inclined to believe presumably that actual ex- 
istence of several non-dipole fields at that time. 
If these assumptions would be accepted, one 
may expect that the geomagnetic poles of the 
same intermediate stage could in more cases 
be found to locate at random in the low lati- 
tudes. 


KOBAYASHI K. (Faculty of Liberal Arts 
and Science, Shinshu University, Matsumoto) 
Geomagnetic field about the time of 
the Plio-Pleistocene boundary (1) Trans- 
itional variation of geomagnetic field since 
Pliocene times hitherto ascertained in Shinshu 
district of Central Japan, are expressed as: 
Ne>Rw—>Rw—Re>N. Where fw indicates the 
intermediate direction. It is impossible, how- 
ever, to place the border line between the 
Pliocene and the Pleistocene at any horizon 
within this Cenozoic section, basing fossil evi- 
dences. (2) As many geologists are acquainted 
with, the demarcation between the Pliocene and 
the Pleistocene was fixed at the meeting of 
International Geological Congress held in 1948 
at London, resulting that the supposed bases 
of the Pleistocene have been improved in many 
parts of the world. (3) On this line of dis- 
cussions, it seems to be needed that a standard 
succession representing geomagnetic variations 
could be established in the region where the 
stratigraphic horizon equivalent to the base of 
the Villafranchian or that being safely correlated 
with the type section in northern Italy could 
be adopted. (4) A geological correlation chart 
was made basing many writings, especially by 
Venzo, Movius, Boule, Haug, Woldstedt and 
others. Moreover the geomagnetic directions 
known in the sense of the normal and the 
reverse were put side by side in a chart to 
make a correlation with the former. The sources 
were found in the writings by Hospers, Einars- 
son in Iceland, Roche in France, Nagata et al., 
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Shimizu, Kawai and Momose in Japan. (5) The 
knowledges afforded by Roche who made his 
works upon the rock sample taken from the 
basin Puy-en-Velay in France where the affairs 
being favourable to correlate the beds with the 
type section of northern Italy, are in their 
reliability the most usable. (6) As far as I 
know, the geomagnetic directions during the 
Astian and the very base of the Villafranchian 
are needed for future measurements. (7) In 
our area of Central Japan, the method of corre- 
lation by means of rock magnetism have 
supplied a lots of valuable knowledges to geology 
of volcanoes perhaps of later Pliocene times. 
Although for the purpose to make use of geo- 
magnetic variation as a standard of demarcation 
of the Plio-Peistocene boundary, the matters 
seem to await further studies, the method will 
offer a global scope of correlation in which we 
might not only settle this tough problem but 
also examine the bases and validities of geo- 
logical correlations. 


NAGATA T., AKIMOTO S., SHIMIZU Y., 


KOBAYASHI K., SYONO Y. — (Geophysical 
Institute, Tokyo University) and KUNO H. 
(Geological Institute, Tokyo University) 


Palaeomagnetic Study on Cretaceous and 
Jurassic Rocks in Kitakami District 
Results of palaeomagnetic studies on Cretaceous 
and Jurassic igneous and sedimentary rocks in 
Kitakami district, North-eastern Japan, are 
summarized in the table, where specimens No. 


1801, 2001, 2105, 2106 are assumed to have 
reverse NRM. Several of the measurements 
were made by high-sensitive astatic magnet- 
ometer recently constructed at the Geophysical 
Observatory in Kakioka. As seen in the table, 
the positions of virtual pole at lower Cretaceous 
and Jurassic ages are situated in southern 
hemisphere.. Especially it is remarkable that 
those from Jurassic samples are in fairly high 
latitude in southern hemisphere. 


NAGATA T. and YAMA-AI M. (Geo- 
physical Institute, Tokyo University) Palaeo- 
magnetic studies on rocks from the coast 
of Liitzow-Holm Bay Results of palaeo- 
magnetic studies on Cambrian gneiss of Ongul 
Islands have been already reported by Nagata 
and Shimizu. During the wintering period of 
JARE III, rock samples were collected at various 
localities along the coast of Liitzow-Holm Bay 
as well as at East and West Ongul Islands. 
Ferromagnetic minerals carrying NRM contained 
in the rocks are pyrrhotite and magnetite with 
small amount of hematite, while those of the 
rocks used in the previous report are mostly 
hematite. It is shown from the measurement 
of NRM that the present rock-samples can be 
classified into two groups, namely, normally 
magnetized one and reversely magnetized one. 
Average values of direction of NRM of these 
two groups of Ongul Islands and Coast of 
Liitzow-Holm Bay are —21°D, —69°I; —10°D, 
49°I and —39°D, —83°I; —19°D, 60°I respec- 

— 


Table 
Specimen Locality Age eee D* jx | Pole Error |No. of Intensity 
No. Position | Angle Specimens] range si 
gr. 
lower 
1801 |Ohunato | Creta- | tuff N 54°W]  0°|152°W 27°S| 6° 290 0-2~ 2.45 
< ceous en 
basalt e : : 
1802 |Kesennuma) | asalt | |N179°W| 46°|141°E 24°S| 10° 16... ae ae 
ee andesiti : ~ i ; 
2001 | Osima 1 vHlava, 61° W) “to* 142° Wr 25ee Gs eee ces 
2002 |Osima Ha | N 74°W|— 6°|142°W 10°S| 9° | 15 3-14~ 29.0 
x10 
2003 | Osima II black ° el yane 
sima Tb pane uie|N 83°W|—16°|140°W 11°S| 8° a a a 
\ 
2105 ‘| Simi ic| dyke ° : 
Simizuhama| Jurassic shale N 33°W| 41°|147°W 57°S| 14° 13 ~10-6 
Cape of sandst 
2106 |Simizuhama| " |shale  |N 26°W| 34°|155°W 59°S| 12° 9 ~10-8 


* after tilt correction 
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tively before tilt correction. As the direction of 
the geomagnetic field at locality concerned here 
is about —43° in D and about —66° in I, it 
may be said that the mean direction of mag- 
netization of the normal Samples is in rough 
agreement with that of the present geomagnetic 
field. Those normal samples are magnetically 
weak and rather unstable. Thus, it may be 
reasonably concluded that NRM of the normal 
Samples is due to accumulation of IRM vis- 
cously caused by the present field, and con- 
sequently that these must be eliminated in the 


present palaeomagnetic studies. North pole 
positions corresponding to the magnetic orien- 
tation thus obtained in the previous and present 
Studies are as follows; 19°N, 163°W (East Ongul 
Island by N. and S.), 9°N, 148°W (East and 
West Ongul Islands), 21°N, 156°W (East Coast 
of Liitzow-Holm Bay). Such agreement of the 
results from three independent sets of sample 
collections at different localities seems to indi- 
cate that the conclusion reported in the previous 
report is fairly reliable. 
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